On the Self-Affine Roughness of a Crack Front in Heterogeneous Media 
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The long-ranged elastic model, which is believed to describe the evolution of a self-affine rough 
crack-front, is analyzed to linear and non-linear orders. It is shown that the nonlinear terms, while 
important in changing the front dynamics, are not changing the scaling exponent which characterizes 
the roughness of the front. The scaling exponent thus predicted by the model is much smaller than 
the one observed experimentally. The inevitable conclusion is that the gap between the results of 
experiments and the model that is supposed to describe them is too large, and some new physics 
has to be invoked for another model. 



The self-affine roughness of a crack-front propagating 
under a tensile load in a randomly heterogeneous system 
is a well studied issue, both experimentally and theoret- 
ically. Experimentally one measures the position h(x, t) 
of the crack front, where h is the position of the front as a 
function of the span- wise coordinate x at time t, and finds 
that this is a self-affine function whose roughness is char- 
acterized by a scaling exponent C (defined below in Eq. 
(|2])) in the range of 0.5-0.65 [j], 0]. Theoretically there 
appears to be a consensus that the appropriate model for 
such dynamical roughening is a long-ranged elastic string 
close to its depinning threshold. This model is defined 
by the equation of motion for a front h(x,t), which is 
allowed to move only forward due to the irreversibility of 
the fracture process Q 
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Here and bellow the integral is meant in the Cauchy Prin- 
cipal Value sense. The RHS of Eq. |T]) is the difference 
between the local driving force (below referred to as G, 
related physically to the energy release rate driving the 
crack [J]), and T(x,h) which is a random quenched noise 
(representing the random material fracture energy [ij]). 
G(°) is the control parameter that represents the energy 
release rate of a straight front. The integral term stands 
for the long ranged restoring forces stemming from bulk 
elastic degrees of freedom. The correspondence between 
the theoretical model and the experimental findings re- 
mained however unclear, since the best numerical stud- 
ies of the resulting self-affine graph h(x) of this model 
came up with a roughness exponent £ = 0.388 ± 0.002 
Q , clearly outside the range of error of the experimental 
measurements. This apparent difficulty led to a number 
of interesting studies, insisting that the model is basically 
right, and that the result concerning the scaling exponent 
is not final. Thus, for example, in the authors ana- 
lyzed Eq. ([1]) using a functional renormalization group. 
They have calculated the scaling exponent £ to one- and 
two-loop orders in e, where e = 2 — d. To one-loop or- 
der the result is ( = e/3, predicting £ = 1/3 at d = 1, 
deviating considerably from the best numerical estimate. 
To two-loop order the prediction increases to about 0.466, 



leading to a statement that the model probably describes 
properly the experimental findings. Unfortunately, it is 
well known that the e expansion is often an asymptotic 
series 0] , sometime providing a better estimate of the ex- 
ponents at first order than at second order. For all that 
one knows the third loop-order may bring the exponent 
down, maybe even below the first loop-order. Another 
approach was advocated in Refs. [1, who proposed 
that the discrepancy between model and experiment may 
be assigned to the existence of nonlinear contributions to 
Eq. In I10| the authors derived, in agreement with 
the results of II ill . that to second order in nonlinearity 
Eq. (JU reads 0|: 
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= J_ ff [h(x\t) - h(x,t)][h(x",t) - h(x',t)} „ 
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where the prime denotes a derivative with respect to x. 
Note that the coefficients are all determined by elasticity 
theory and are not free. Using a method proposed by 
Schwartz and Edwards fl3| . it was concluded that the 
nonlinear term affects the scaling exponent dramatically, 
stating that £ > 0.5 @. On the other hand in Q a similar 
nonlinear equation was analyzed in the framework of one- 
loop renormalization group, yielding £ 0.45. Clearly, 
the situation warrants some further scrutiny. 

In this Letter we present careful numerical measure- 
ments of the scaling exponent of the present model to 
first and second order in h(x,t). To this aim we simulate 
the dynamical model |T]) with and without the nonlinear 
terms in Eq. @. Our final conclusion is that although 
the second order terms perturb the solution h(x, t) signif- 
icantly, they are actually irrelevant for the scaling expo- 
nent, that appears unchanged with or without the non- 
linear terms. The uneasy conclusion of this analysis is 
that the model itself may not describe the experiment 
correctly; a discussion of this conclusion is offered at the 
end of this Letter. 

To numerically simulate the model we discretize the 
spatial variable x, and swap temporal changes with dis- 
cretized steps in the variable h(x). Choosing the ba- 
sic unit of length to be in the x direction, we present 
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below simulations for x £ [1,L] with L = 2™ in the 
range [2048,16384]. We used periodic boundary condi- 
tions. The discretization of h(x) is chosen in units of 1/7; 
this seems arbitrary, but since in the depinning limit the 
velocity is irrelevant, this discretization should not af- 
fect the scaling exponents. At step zero the interface is 
prepared with a random jitter to avoid spurious lattice 
artifacts. The random quenched noise T(x,h), is picked 
at each lattice point from a uniform distribution in the 
interval [0,1.5] [lij]. Following [l5| we simulate the de- 
pinning limit by increasing G^ incrementally from zero, 
such that the local driving force at the least pinned point 
overcomes the local fracture energy T. This local depin- 
ning may trigger additional motion until the interface 
stops, at which moment G^ is increased further until 
the next weakly pinned point gives in. Measurements of 
the roughness were taken when all the points x moved at 
least one step after the last increment in G^°\ 

For the present calculation we employed the rms defi- 
nition of the roughness, i.e. 
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FIG. 1: Demonstration of data collapse for the linear model, 
plotting log 10 O(^,L)L" c ) vs. \og w (£/L) with L = 2", n=12, 
13 and 14, using £ = 0.362. The results are obtained by 
averaging over 100 realizations. 
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We keep the implicit L dependence in this quantity since 
it turns out that the roughness exponent is a slowly con- 
vergent quantity as a function of L. Indeed, it was con- 
vincingly demonstrated in Q that the numerical value of 
roughness exponent as measured using the linear model 
(fTJ converges in the relation 



w(£,L) ~^ 



(4) 



only when L is of the order of 10 6 . Not having simulations 
of this order we resort to finite size scaling which was 
demonstrated [l6| to yield reliable exponents also with 
smaller values of L. The essence of this method is the 
finite-size scaling assumption written as 



w(i,L)=L<f{l/L) 



(5) 



for L very large there is a range of values of £ where 
f(l/L) ~ (l/L) 1 *, coalescing with the simple scaling as- 
sumption For smaller values of L one seeks the best 
value of ( by data collapsing w(£, L)L~^ onto a univer- 
sal function f(£/L). An example of this procedure is 
shown in Fig. [I] in which the first order model had been 
employed and the results were averaged over 100 realiza- 
tions for each L. The data collapse appears satisfactory 
with the choice ( — 0.362. The degree of confidence that 
this method provides can be demonstrated by the inferior 
data collapse obtained for the same data with ( = 0.4 in 
Fig. O Our best estimate of the scaling exponent of the 
model realized to first order is £ = 0.365 ± 0.005. Note 
that this exponent is higher than the value £ ~ 0.35 
obtained using Eq. ((3]) with L — 16384. This is in agree- 
ment with the statements in the literature for the slow- 
ness of convergence of the scaling exponent with L p|. 
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FIG. 2: Demonstration of the failure of data collapse for the 
linear model, plotting log 10 («j(£, L)L~^) vs. log 10 (£/L) with 
L = 2 n , n=12, 13 and 14, using £ = 0.4. The results are 
obtained by averaging over 100 realizations. 



The finite size scaling analysis improves the situation, 
even though our estimate still falls a bit short compared 
to the estimate C = 0.388 for L w 10 6 This difference 
will not pose a difficulty in assessing the importance of 
the nonlinear term. 

Adding the nonlinear terms, one should first ascertain 
that they make a significant change in the front dynam- 
ics. This is demonstrated in Fig. [3l which compares, for 
the same initial interface and the same quenched noise 
r(x, h{x)) 1 the realized interfaces, once with only a lin- 
ear term and once with the full second-order nonlinear 
contributions. It is obvious that the nonlinear terms are 
not negligible, they inflict major changes on the actual 
graph. The seasoned reader will notice however that the 
scaling exponent is hardly changed. This eye-ball con- 
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FIG. 3: Comparison of the fronts obtained with the linear 
and the nonlinear model, where the initial condition and the 
quenched noise are all the same. Note the huge difference in 
scales between the abscissa and the ordinate. 
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FIG. 4: Demonstration of data collapse for the nonlinear 
model, plotting log 10 (w(£, vs. log 10 (^/L) with L = 

2 n , n=12, 13 and 14, using £ = 0.362. The results are ob- 
tained by averaging over 100 realizations. 

elusion is fully supported by the finite-size scaling anal- 
ysis which is presented in Fig. [4] The quality of the 
data-collapse is essentially identical in Figs. [1] and 2] its- 
ing the same scaling exponent in both cases. Note that 
the amplitude of the overall front roughness is reduced 
in comparison to the linear model; the nonlinearity in- 
creases the stiffness of the front. The exponent remains 
however invariant. We thus conclude that the numerical 



evidence presented here does not support the theoretical 
propositions of [j| Q . 

In light of these results we propose that the relation be- 
tween the experiments and the model must be reassessed. 
One could argue that our fronts are not large enough to 
asymptote to a "correct" scaling behavior. To such a 
claim one must answer that the typical experiments do 
not have more scales than our simulation. For example 
in the fracture experiments of [l[ the randomness scale 
(also known as the correlation length) is determined by 
the size of sand particles that blast the interface between 
two slabs of material that are then glued together. This 
typical scale, which determines the scale of the fracture 
energy T in the present model, is of the order of 50^m. 
Crack front segment up to 50 mm were analyzed, giving 
at most three orders of scales in theory, but in practice 
the self affine scaling was observed in a range that spans 
about two orders of magnitude. The measured roughness 
exponent is significantly larger than the values discussed 
above, even in this limited range of length scales. It is 
therefore entirely reasonable, in our opinion, to restrict 
the theoretical analysis of any given model to about the 
same range of length scales or slightly more, as is done 
above. Theories invoking asymptotically large system 
sizes may be interesting, but hardly relevant for such ex- 
periments. 

Accordingly we may ask what is missing in the relation 
between theory and experiment. One thing that may be 
suspicious is the assumption at the background of the 
derivation of the models (fT|) and ©, i.e. that elasticity 
theory is entirely sufficient to describe the crack front 
dynamics. Since elasticity theory predicts the divergence 
of stress at the crack front [i[ , realistic materials will al- 
most surely yield plastically or develop additional local 
damage. Such a change in material properties, precisely 
where the dynamics is taking place, may very well change 
the nature of the long-range interactions of the bulk de- 
grees of freedom. How to renormalize such long-range 
interactions when plasticity or other modes of damage 
are at play is not known at this point in time. We stress, 
however, that the gap between the model results and 
the experimental results indicate that such novel thinking 
about the theoretical fundamentals may be unavoidable. 
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